Solar Physics 
DOI: 10.1007/' 



On the Relation between Solar Activity and 
Clear-Sky Terrestrial Irradiance 

G. Feulner 

© Springer •••• 

Abstract The Mauna Loa Observatory record of direct-beam solar irradiance 
measurements for the years 1958-2010 is analysed to investigate the variation 
of clear-sky terrestrial insolation with solar activity over more than four solar 
cycles. The raw irradiance data exhibit a marked seasonal cycle, extended periods 
of lower irradiance due to emissions of volcanic aerosols, and a long-term decrease 
in atmospheric transmission independent of solar activity. After correcting for 
these effects, it is found that clear-sky terrestrial irradiance typically varies by 
~ 0.2 ± 0.1% over the course of the solar cycle, a change of the same order of 
magnitude as the variations of the total solar irradiance above the atmosphere. 
An investigation of changes in the clear-sky atmospheric transmission fails to hnd 
a significant trend with sunspot number. Hence there is no evidence for a yet 
unknown effect amplifying variations of clear-sky irradiance with solar activity. 

Keywords: Atmospheric extinction; Integrated Sun observations; Solar cycle, 
observations; Sunspots, statistics 



1. Introduction 

Long-term records of clear-sky terrestrial so lar irradiance are very use ful to study 



changes in atmospheric transmission fe.q.. lHovt and Frohlichl . 119831 ). It is more 
difficult to derive changes in insolation due to solar variability from these data 
because the variations related to solar activity are much smaller than the changes 
in atmospheric transmission of different origin and the systematic and random 
errors in the data. Nevertheless, such an approach has been tried many times in 
the past. 

The Smithsonian Astrophysical Observatory (SAO) observations from the 
first half of the 20th century are an important historical example for terrestrial 

i rradiance measurements with the aim to detect solar variability from the ground 

(| Abbot. Aldrich. and Fowld . Il932t I Abbot. Aldrich. and Hooverl Il942l : lAldrich and Hooverl .| 



G. Feulner 

Earth System Analysis, Potsdam Institute for Climate 
Impact Research (PIK), P.O. Box 6012 03, 
D-14412 Potsdam, Germany 
e-mail: feulner@pik-potsdam.de 



SOLA: maunaloa.tex; 10 October 2012; 1:32; p. 1 



Fculncr 



1954I ). The SAO observing campaigns were designed to derive a time series of 
the solar constan t or total so lar irradiance (TSI) above the atmosph ere, but 
calibration issues (lAllenl J 19581) , biases due to strong seasonal variations (jFeulnerl 



l201llaL and chan ges in atmospheric transmission independent of solar a ctivity 
(IAngstr5ml 1970|) - mostly relate d to volcanic aerosols and local pollution ( Hovt 



19791 : Roosen and Angionei . 1984 ) - can pose serious difficulties for investigations 



of irradiance changes related to solar activity if not properly accounted for 



De spite t hese problems, the SAO data continue to attract attention. IWeber 



( 2010l . 2011 ) correlated the ground-based irradiance measurements from the SAO 



data with sunspot numbers and claimed to have found a strong variation of ter- 
restrial insolation with solar activity, with differences of the order of 1% between 
solar maxima and minima (i.e., one order of magnitude larger than the v ariations 
in the TSI on top of the atmosphere) . These findings were criticised by iFeulner 



(2011a, b) who could show that they were due to seasonal bias and the effects of 
volcanic aerosols and local pollution. The seasonal bias arises from the fact that, 
by coincidence, days with large sunspot numbers in the SAO data predominantly 
lie in months with seasonally high atmospheric transmission, making irradiance 
values at times of high solar activity appear higher. The second bias is introduced 
because two out of three solar minima in the SAO data covering solar cycles 16, 
17, and 18 are affected by reduced atmospheric transmission due to volcanic 
aerosols and local pollution. This effect makes solar irradiance on the ground 
during those two minima appear lower, thus again suggesting a stronger variation 
of terrestrial irradiance with solar activity. After correcting for the seasonal cycle 
and excluding periods of time affected by aerosols from volcanic eruptions and 
local pollution, the variations of terrestri al insolation arc of the same order of 
magnitude as the TSI variations dFeulnerl. 2011a. b). 



More recently. Irlempelmann and Weberl (|2012[ ) used terrestrial insolation data 
taken at Mauna Loa Observatory (MLO) and reported variations of the solar 
irradiance on the ground a factor of 10 larger than above the atmosphere. Here 
an improved analysis of these data is presented to independently quantify how 
much clear-sky terrestrial irradiance varies between solar maxima and minima. 

This paper is organised as follows. Section [2] describes the Mauna Loa terres- 
trial insolation data and corrections of the effects of the seasonal cycle, aerosols 
from volcanic eruptions, and a linear long-term trend independent of solar ac- 
tivity. In Section [3] the correlation of the thus corrected irradiance with sunspot 
number is presented. Section U analyses the atmospheric transmission factor and 
its changes with sunspot number, before the results arc discussed in the context 
of previous work in Section [5] Finally, Section [B] summarises and concludes this 
paper. 



2. Data Analysis 

2.1. A First Look at the Raw Irradiance Data 

For the analysis in this paper, the ground-based direct-beam solar irradiance 
data taken since 1958 at Mauna Loa Observatory (MLO, latitude 19.533° N, 
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longitude 155.578° W, elevation 3400 m) are used fellis and Pueschell 11971 



Dutton. Deluisi. and Austring. 1985; Dutton et all 1994; D utton and Bodhaine 



20011 ). These measurements have been obtained at local noon as well as at 
airmass values of 2, 3, 4, and 5 (both in the morning and in the afternoon). In 
the following, results for the morning irradiancc data at airmass 2 are presented. 
Airmass 2 is chosen since smaller airmass values correspond to smaller zenith 
angles and thus shorter paths through the atmosphere. Furthermore, as in most 
previous studies based on these data, the morning measurements are preferred 
over the afternoon data or a combined record of the two due to the absence 
of local influences on th e atmosphe r ic tra nsmission because of downslope winds 
in the morning hours ( Mendoncal . 1969b . Indeed, a comparison of irradiance 
morning and afternoon data at airmass two for days where both measurements 
are available shows that morning irradiance values are larger on average and 
thus apparently less affected by local disturbances in atmospheric transmission 
(see Figured]). 
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Figure 1. Difference of morning and afternoon irradiance measurements at MLO and airmass 
2 for days where both are available (red histogram) , showing that morning values are on average 
larger than afternoon values. The blue curve shows a Gaussian function fit to the data with a 
peak at 12.9 W m" 2 and a width cr = 17.29 W m~ 2 . 



As a consistency check, the full analysis has been carried out for a combined 
record of morning and afternoon data and for the other airmass values as well, 
however, yielding very similar results which are not shown in this paper. 

The raw irradiance data at morning airmass 2, for historic reasons given in 
units of Langley min -1 or cal cm~ 2 min , arc shown in Figures [5^, and fjZb- 
In this figure, they are compared to the globally averaged op tical depth trrq o f 
volcanic aerosols at a wavelengt h of A = 550 nm as com piled bv lSato et all ( 1993f) 
for the years 1959-1993 and bv lSolomon et all (|201ll ) for 1994-2010 (FigureG_p) 
and to daily international (also called Wolf or Zurich) sunspot numberaj R 
(Figured;). 



Available at http://sidc.oma.be/DATA/dayssn_import.dat, date of access: 22 December 2011. 



SOLA: maunaloa.tex; 10 October 2012; 1:32; p. 3 



Fculncr 



2.0 
1.5 

'E 1.0 

CM 

| 0.5 
"S3 

o, o.o 

o 1-8 
g 1.6 

"3 1-4 
J= 1.2 
1.0 
0.8 

o 

* 0.15 
_c 

Q. 

_2> 0.10 
Q 



-(b) 

| iflvvi 

— ■"' ? 



g 0.05 h 



300 - 



o 

w 100 



(a) 



A 



ttttttt- , - 7 i rr . -TTT-T | -" v , T~rr i - , •••}....■, ■ 



i 



i ^ Agung Fernandina Kamchatka El Chichon 



I 



Pinatubo 




! (d) 



19 



20 



»U i, 



21 



22 



23 



24 




1960 1970 1980 1990 

Year 



2000 



2010 



Figure 2. The MLO terrestrial irradiance data for the time period 1958-2010 in the context 
of solar activity and volcanic eruptions during this time interval, (a) Raw irradiance measure- 
ments for morning measurements at airmass 2. (b) Same as (a), but zooming int o the data, 
(c) G lobal average optical depth rgso o f volcanic aero s ols at A = 550 nm from ISato et al\ 
(1993) until the end of 1993 and from ISolomon et ali J2011P| since the beginning of 1994. 
Prominent volcanic eruptions are marked, (c) Daily international sunspot numbers R during 
this time interval, spanning a large part of solar cycle 19, cycles 20, 21, 22, and 23 as well as 
the beginning of the current cycle 24. 



Several interesting observations can be made in the raw data which are 
relevant for a proper analysis of correlations with solar activity: 



i) The irradiance data exhibit a pronounced seasonal cycle. As shown in lFeulner 
(2011a), not correcting for seasonal variations can affect correlations with 
solar activity if different phases of solar activity are not evenly distributed 
over the seasons in the data. 

ii) There is a cluster of irradiance data around 0.5 cal cm -2 min" 1 in the years 
1975-1976. As one can clearly see the seasonal cycle and as these values are 
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offset by 1 from the majority of the data; these are most likely due to typos 
in the record where values have a leading digit of instead of 1. These typos 
can either be corrected or masked out; the latter option has been chosen in 
this paper. 

Hi) The dimming effects of volcanic aerosols in the atmosphere can be seen di- 
rectly in the data. These are most pronounced for the largest eruptions, i.e., 
El Chichon in 1982 and Pinatubo in 1991, of course, but irradiance values 
are also lower during volcanic episodes in the 1960s and 70s. This can be best 
seen in the expanded view presented in Figure [5Jd. 

iv) Looking at periods only weakly affected by volcanic aerosols over the entire 
record, one can see a slight decrease of te rrestrial irradiance with time pre- 
viously described in lSolomon et al. (l201lh . It will be shown below that this 



long-term trend is well approximated by a linear fit and not connected to 
changes in sunspot number. 

All these issues have to be taken care of before correlating the MLO ter- 
restrial irradiance data with sunspot number in order to ensure an unbiased 
estimate of the influence of solar activity on clear-sky irradiance on the ground. 
In addition, the MLO data are also affected by changes in instrumentation 
and offsets in calibra t ion w hich can also influence any trend analysis (e.g., 



Dutton and Bodhaind . l2001r ) . These effects can be el iminated, however, by di 



rectly investigating atmospheric transmission factors (jEllis and Pueschell 119711 ) 
as shown in Section 0] 



2.2. Corrections Applied to the Irradiance Data 



As a first step, the raw MLO irradiance data for morning airmass 2 as shown in 
Figures [2^, and[2jD are converted to SI units (W m -2 ) by multiplying by a factor 
of 697.4. These converted irradiance data are shown both as a time series and as 
a scatter plot vs. sunspot number in Figure [3^,. Furthermore, for each day with 
observations the MLO irradiance data are correlated with daily international 
sunspot numbers and interpolated volcanic aerosol optical depth data from the 
same day. 

In the discussion above, three sources of systematic bias of the irradiance data 
were identified: the seasonal cycle, attenuation by volcanic aerosols, and a pos- 
sible long-term trend. These effects will be discussed and, if possible, corrected 
in the following. The individual corrections and their effect on the time series 
as well as the correlation between irradiancc and sunspot number are shown in 
Figure 



2.2.1. Seasonal Cycle 



The pronounced seasonal cycle is one of the most obvious features of the raw 
irradiance time series shown in Figure I3J}. One way to correct for this effect is to 
take data from several years only little affected by volcanic aero sols and combine 
them to compute median monthly values for the irradiance (jFeulnerl . 2011a). 
Daily values of this average seasonal irradiance cycle can then be computed 
from a cubic-spline fit through the data with periodic boundary conditions. 
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Figure 3. Corrections applied to the MLO terrestrial irradiancc data at airmass 2. The 
left-hand panels show the irradiance time series, the right-hand panels scatter plots of irradi- 
ance vs. sunspot number, (a) Raw irradiance data as shown in Figures [2ji and[2p. The grey 
shading indicates the portion of the data used to derive the average seasonal cycle in terrestrial 
irradiance. (b) The raw data with the seasonal cycle subtracted, (c) The data corrected for 
both seasonal cycle and for the effects of volcanic aerosols using the globally averaged optical 
depth shown in Figure [2p. The two-year time periods with residuals after the El Chichon and 
Pinatubo eruptions excluded in the analysis below are indicated by the grey shading. The red 
line shows the linear long-term trend derived from a fit to the data, (d) Same as (c), but with 
the long-term trend subtracted. 



For the MLO data the years from 2000 to 2009 have been chosen to compute 
the average seasonal cycle (indicated by the grey shading in Figure [3^) as this 
period of time is only weakly affected by volcanic aerosols. Furthermore, a linear 
trend with time has been subtracted before computing the monthly median irra- 
diancc values to remove any long-term changes in the irradiance baseline (be it 
from solar activity changes, trends in volcanic or other aerosols, or instrumental 
drift). The resulting cycle of seasonal irradiance anomalies is shown in Figure [4l 

The irradiance data with the seasonal cycle anomaly removed are shown in 
Figure The irradiance variations through the year are clearly dominated by 
the changes in the distance to the Sun on Earth's elliptical orbit. In addition 
there are minor, but statistically insignificant, seasonal anomalies, in particular 
in spring and summer. 

Note that there is quite a bit of intcrannual variation in the seasonal irradiance 
cycle; therefore, a perfect correction of the effect cannot be expected. Comparing 
the irradiancc time series with the seasonal cycle subtracted shown in Figure [3Jd 
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Figure 4. Average seasonal cycle in the terrestrial irradiance at MLO and airmass 2 computed 
from the irradiance data in the years 2000-2009 (small black squares) and after removing a 
linear long-term trend. The filled red circles are the monthly medians while the red line shows 
a fit of cubic splines with periodic boundary conditions through the monthly medians. The 
blue line indicates the expected seasonal irradiancc variations due to Earth's orbital motion. 



to the raw irradiance data shows, however, that the seasonal cycle is effectively 
removed. In particular, the sequence in the irradiance vs. sunspot number scatter 
plot shown in the right-hand panel of Figure exhibits less scatter and visually 
appears considerably better defined. 

2.2.2. Effects of Volcanic Aerosols 

The second atmospheric effect which has to be taken into account when in- 
vestigating the correlation between solar activity and terrestrial insolation are 
the attenuating effects of volcanic aerosols. One way of dealing with episodes 
of volcanic eruptions is to exclude these periods fro m further analysis as has 
been done for the SAO terrestrial irradiance data in Feulnerl (2011a). This is 



not a viable option for the frequent volcanic eruptions affecting the MLO data, 
however. The expected irradiance variations due to solar activity are below 1% 
or AI < 10 W m~ 2 corresponding to an attenuation by volcanic aerosols with an 
optical depth of t < 0.01. Such a cut in aerosol optical depth would mean that 
more than half of the MLO irradiancc record would have to be excluded from the 
analysis, and even then the remaining record would suffer from inhomogencous 
attenuation by volcanic aerosols resulting in irradiance variations at least of the 
order of the changes expected from solar variability. 

A preferable option would be to correct for the effects of volcanic eruptions 
using an independent da taset of the optica l depth of volcanic aerosols. For the 



analysis in this paper, the lSato et all ( 19931) compilation of global average optical 



depth T550 of volcanic aerosol s at wavelength A = 5 50 nm is used for the years 



1958-1993, extended with the lSolomon et all (| 20111) data until March 2010 and 



with a constant value T550 = 0.007 until the end of 2010. The irradiance I s , 
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corrected for the seasonal cycle and the effects of volcanic aerosols is then cal- 
culated from the irradiance I s , which is corrected for the seasonal cycle only, 
using I sv = I s exp(X cx T550). In this formula, X denotes the airmass, and cx is 
a correction factor converting the optical depth at 550 nm to the optical depth 
integrated over the whole visible solar spectrum. The correction factors cx for 
the four airmass values are determined empirically, finding C2 = 0.50 ± 0.01, 
c 3 = 0.45 ± 0.01, c 4 = 0.42 ± 0.01 and c 5 = 0.40 ± 0.01. 

The time series and scatter plot for the irradiance / sv corrected for seasons 
and volcanic aerosols is shown in Figure Overall the correction seems to work 
well for the late stages of volcanic eruptions. There are, however, considerable 
residuals during times of rapidly changing volcanic aerosol load, in particular 
for the Pinatubo and El Chichon eruptions. This is not surprising since the 
global average aerosol optical depth may not be fully representative for the local 
aerosol load abov e Mauna Lo a , in pa rticular in the period following the eruptions. 



Furthermore, the ISato et all (|19931 ) compilation provides monthly values which 
are not able to trace the fast changes visible in the daily irradiance record. Two- 
year periods after the four strongest eruptions (Agung, Fernadina, El Chichon, 
and Pinatubo) will be masked out for the analysis of the correlation between 
terrestrial irradiance and sunspot number presented below. Note that there is 
also a dip in the thus corrected irradiance around the year 1967. corres ponding 



to a decrease in volcanic aerosol optical depth in the ISato et al\ (|1993f ) record. 
As the origin of this period of lower irradiance remains unclear, it will not be 
excluded from the analysis. 

2.2.3. Long-Term Trend 

Looking at the scatter plot of irradiance values corrected for the seasonal cycle 
and volcanic aerosols with sunspot number shown in the right-hand panel of 
Figure [3J: one can indeed see a marked increase of terrestrial solar irradiance with 
sunspot number. At least part of this trend is clearly driven by what appears to 
be a separate sequence of measurements at high irradiance (I ~ 1100 W m~ 2 ) 
at relatively high sunspot numbers (100 < R < 300). 

A comparison with the irradiance time series shows that the vast majority of 
these measurements are from the very first years of the record (years 1958-1961) 
which, according to Figure [21 arc unaffected by volcanic aerosols and coincide 
with the particularly strong maximum of solar cycle 19. 

These observations are important since the irradiance data shown in Figure [3fc 
exhibit a slowly decreasing trend over time. This general decrease in irradiance 
can be clearly seen during the last decade, for example, or by comparing the 
beginning of the record and the end. This trend in atmospheric transmission 
has been not ed before and a t tribut ed to changes in the atmospheric background 



aerosol load (jSolomon et al\ . 120111 ) 



It will be shown below that this long-term trend can be detected in the 
data spanning almost five solar cycles. Note, for example, that the decrease 
can be seen over the entire solar cycle 23 with a constant slope despite large 
variations in sunspot number. The trend can thus be regarded as independent 
of solar activity and has to be subtracted to ensure an accurate analysis of the 
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correlation between solar activity and terrestrial irradiance. To this end a linear 
trend It (t) = Iq + ai t is approximated to the entire record of irradiance data 
corrected for the seasonal cycle and volcanic aerosols, but excluding the two- 
year time intervals after the El Chichon and Pinatubo eruptions. Furthermore, 
the data used to fit the linear trend are restricted to the range < R < 200 
to ensure that the fit of the long-term trend is not driven by changes in solar 
activity between the strong solar maximum of cycle 19 and the considerably 
less pronounced maximum of cycle 23. Fitting a linear trend using these data 
yields a slope of qj = — 0.38 ± 0.02 W m -2 yr _1 . This is very similar to the 



trend found in ISolomon et all (|201lh . As will be shown below, a multivariate 
regression has been performed on / sv to simultaneously analyse its linear trends 
with sunspot number and time as an additional test whether the trend in time 
reported here is indeed independent of solar activity, finding a very similar value 
for the long-term decrease in terrestrial irradiance with time. 

Thus there is indeed a significant long-term decrease in terrestrial irradiance 
which is independent of solar activity and therefore has to be subtracted from 
the data before analysing any correlation of irradiance with sunspot number. 
The resulting irradiance time series of I sv t(t) = Isv(t) — h(t) and its scatter 
plot are shown in Figure [3ji. It should be noted that the conspicuous second 
sequence at large irradiance values in the scatter plot has now disappeared and 
that the overall correlation shows a considerably smaller slope than before the 
subtraction of the linear long-term trend. 



3. Correlation of Terrestrial Irradiance with Sunspot Number 

After applying the corrections for the seasonal cycle, for the attenuating effects 
of volcanic aerosols, and for the linear long-term trend, the MLO terrestrial 
irradiance data at morning airmass 2 can now be correlated with the sunspot 
number to investigate changes of terrestrial irradiance with solar activity. 

The corrected irradiance data J sv t as a function of sunspot number R are 
presented in Figure Fitting a line to this correlations yields a slope 



disvt/di? = +(0.019 ± 0.007) Wm" 



(1) 



The error has been computed from 10 000 bootstrapping simulations where the 
sample is first duplicated and added to the original sample, then modified by 
random errors, before half of the enlar ged sample is randomly selected to perform 
the linear regression (see also lFeulnerl . 2011a). A ra ndom error of 0.5% is assumed 
for the individual MLO irradiance measurements (jDutton and Bodhaind . 120011 ) 
during these bootstrapping simulations. 

In the analysis so far, the long-term trend with time t of the irradiance 
J 8V corrected for the seasonal cycles and the effects of volcanic aerosols has 
been removed before correlating with sunspot number R. While this is highly 
illustrative, it is statistically more appropriate to perform a simultaneous mul- 
tivariate linear regression of / sv with respect to trends with t and R. Thus 
one can also independently test whether the linear trend in time reported in 
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Figure 5. Scatter plot of MLO terrestrial irradiance vs. sunspot number (black squares). The 
irradiance data are corrected for the seasonal cycle, volcanic aerosols, and a long-term decrease 
in atmospheric transmission independent of solar activity. The red line shows a linear fit to 
the correlation. 



Section I2.2I is indeed independent of solar activity as traced by the sunspot 
number R. This test yields a linear trend of terrestrial irradiance with time of 
dlsv/dt = —0.45 ± 0.02 W m~ 2 yr _1 , which is slightly larger than the value 
reported above, but compatible within the errors. The slope of the correlation 
with sunspot number derived in this way is 



dI sv /dR = +(0.015 ± 0.006) Wm" 2 , (2) 

again in excellent agreement with the value from Equation (1). 

Thus a simultaneous regression of a linear long-term trend with time and 
a linear trend with sunspot number to the terrestrial irradiance data not yet 
corrected for the linear trend yields results which are very similar to the ones 
described above. This is important because if the long-term trend in terrestrial 
irradiance had been caused by changes in solar activity, this long-term trend 
would be reflected in the sunspot numbers as well, and the regression would 
have yielded an insignificant trend with time. This clearly shows that the linear 
long-term trend described above is indeed not caused by solar activity. 

As a test whether this result depends on the choice of sunspot number as 
indicator of solar activity, the simulta neous regression has been repeated fo r 
annual means of the open solar flux F$ (ILockwood. Rouillard. and Finchl . l2009l) . 



finding consistent values for the trends, albeit with larger errors (dl sv /dt = 
-0.47 ± 0.15 W m~ 2 yr" 1 and 8I sv /dF s = (1.4 ± 2.6) x 10~ 14 W mT 2 Wb" 1 ). 

In summary, a positive trend of the corrected terrestrial irradiance data with 
sunspot number can be found in the MLO data which is significant on the ~ 2a 
level. The magnitude of these changes of solar irradiance on the ground with 
respect to changes in irradiance on top of the atmosphere and previous studies 
on correlations of terrestrial irradiance with solar activity will be discussed in 
Section [5j 
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4. Correlation of Atmospheric Transmission with Sunspot Number 
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Figure 6. Same as Figure [3] but for the atmospheric transmission factor (ATF) and the 
corrections applied to this quantity. 



The analysis above investigated the correlation of terrestrial insolation with 
sunspot number to investigate the possibility of a so-far unknown amplification 
effect produced by changes in clear-sky atmospheric transmission. Alternatively, 
one can directly study the atmospheric transmission itself. The advantage of 
this approach is that one can make use of the rationing technique which elim- 
i nates both the variation s in TSI and any instrumental calibration differences 
(jEllis and Pueschell . ll97ll ). In this method, atmospheric transmission is analysed 
in terms of the atmospheric transmission factor (ATF) defined as follows: 



where Ix are the morning irradiance values at the different airmass values X of 
2, 3, 4, and 5. The ATF for the Mauna Loa data is shown in Figure [6^,. 

Before correlating the ATF with sunspot number i?, the same corrections 
have to be applied as for the irradiance data (see Section \2.2\i . The seasonal 
cycle in the ATF is computed as for the irradiance, but note that the dominant 
effect caused by the orbital motion of the Earth is already taken care of by 
the rationing technique. The remaining corrections for the effects of volcanic 
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aerosols and the long-term trend (with a best-fitting value of d(ATF)/d£ = 
— (8.2 ± 0.4) x 10 -5 yr" 1 ) also follow the procedure discussed in detail for the 
irradiance. The stepwise corrections and their effect both on the time series and 
for the correlation with sunspot number are illustrated in Figure [5] 

The correlation between the corrected atmospheric transmission factor ATF 
and sunspot number R is shown in Figure [7J A linear fit to this correlation yields 



d(ATF)/di? = -(0.7 ±1.4) x 10" 



(4) 



As for the irradiance data, an improved analysis is performed in which both 
the long-term trend of the ATF with time and its correlation with sunspot 
number are fit simultaneously. This analysis yields a value of d(ATF)/dt = 
— (10.5 ± 0.4) x 10~ 5 yr" 1 for the long-term trend, in good agreement with 
the value derived above. The fact that the simultaneous linear regression of 
the ATF finds a significant long-term trend with time again confirms that this 
long-term trend is independent of solar activity. The linear correlation of the 
ATF with sunspot number has a slope of <9(ATF)/<9i? = -(2.0 ± 1.3) x 10~ 6 , 
again in agreement with Equation (4) . Thus no significant change of atmospheric 
transmission with sunspot number can be detected in the Mauna Loa data. 




350 



Sunspot Number R 



Figure 7. Scatter plot of MLO atmospheric transmission factor (ATF) vs. sunspot number 
(black squares). The ATF data are corrected for the seasonal cycle, volcanic aerosols, and 
a long-term decrease in atmospheric transmission independent of solar activity. The red line 
shows a linear fit to the correlation. 



To test for a possible dependence of this result on the choice of solar-activity 
indicator, the bivariate regression has been repeated replacing the sunspot num- 
ber R by the open solar flux Fg (using annual averages of the data as for 
the analysis of the irradiance data). Again, the trends derived for the open 
solar flux are in good agreement with the trends computed using sunspot num- 
bers, although the errors are larger (d(ATF) / dt = (-10.7 ± 2.7) yr" 1 and 
<9(ATF)/aF s = (-1.6 ±4.8) x 10" 18 Wb" 1 ). 
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5. Discussion 

The analysis in Section [3] shows that the corrected terrestrial insolation / sv 
linearly increases with sunspot number R with a slope 8I sv /dR = +(0.015 ± 
0.006) W m~ 2 . With typical differences in sunspot number of AR « 150 be- 
tween solar maxima and minima, this translates into irradiance variations of 
AI w 2il W m~ 2 or ~ 0.2 ± 0.1% over the solar cycle. This is of the same 
order of magnitude as the variations in total solar irradiance (TSI) observed 
from satellites (e.g., Frohlich and Lean . 2004 ; Gray et all . l201dh which suggests 



that clear-sky terrestrial irradiance variations with solar activity are not strongly 
amplified by any hitherto un known atmospheric feedba c k pro cesses. 



The conflicting results in iHempelmann and Weberl (|2012l ) who find a one 



order of magnitude larger variation in terrestrial irradiance with solar activity 
(dl/dR = (0.10 — 0.18) W m~ 2 depending on the way the effects of volcanic 
eruptions arc accounted for) can be understood from the fact that the authors 
have not corrected the terrestrial irradiance data for the slow long-term decrease 
in atmospheric transmission not related to changes in solar activity. The com- 
bination of large irradiance values at the beginning of the record with large 
sunspot numbers during the strong maximum of solar cycle 19 then results in 
a spuriously large slope of the irradiance vs. sunspot number correlation. In 
ad dition, excluding periods of up t o three years after volcanic eruptions as done 



in IHempelmann and Weberl (|2012f ) is insufficient as attenuation effects at the 
percent level are visible in the data even after this time, resulting in changes of 
at least the order of magnitude of the changes related to solar variability (see 
the discussion in Section |2~2"|) . 

The lack of evidence for any unknown effect amplifying changes of terrestrial 
irradiance with solar activity in the MLO data is confirmed by an analysis of the 
correlation of the atmospheric transmission factor with sunspot number which 
failed to find a significant change of clear-sky atmospheric transmission with 
solar activity. 



6. Conclusions 

In this paper, the Mauna Loa Observatory (MLO) direct-beam solar irradiance 
record has been analysed to investigate the correlation between solar activity 
and terrestrial insolation. After correcting for the seasonal cycle, the attenuating 
effects of volcanic aerosols, and a long-term decrease in atmospheric transmission 
independent of solar activity, clear-sky terrestrial irradiance is found to vary by 
AI w 0.2 ± 0.1% between maxima and minima of the 11-year solar activity 
cycle. These variations arc of the same order of magnitude as the changes in 
total solar irradiance on top of the atmosphere over the solar activity cycle. An 
investigation of the atmospheric transmission shows that there is no significant 
trend of clear-sky atmospheric transmission with sunspot number. Thus there 
is no evidence for an unknown amplification effect of solar activity in terms of 
strong changes in clear-sky atmospheric transmission. 
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